Fit log returns to F-S skew standardized Student-t distribution.
m is the location parameter.
s is the scale parameter.
nu is the estimated degrees of freedom.
xi is the estimated shape parameter.

Log returns data 2011-2023.

For 2011, medium risk data is used in the high risk data set, as no high risk fund data is available prior to 2012.
vmrl is a long version of Velliv medium risk data, from 2007 to 2023. For 2007 to 2011 (both included) no high risk data is available.

Summary of gross returns

##       vmr             pmr             mmr             vhr       
##  Min.   :0.868   Min.   :0.904   Min.   :0.988   Min.   :0.849  
##  1st Qu.:1.044   1st Qu.:1.042   1st Qu.:1.013   1st Qu.:1.039  
##  Median :1.097   Median :1.084   Median :1.085   Median :1.099  
##  Mean   :1.070   Mean   :1.065   Mean   :1.066   Mean   :1.085  
##  3rd Qu.:1.136   3rd Qu.:1.107   3rd Qu.:1.101   3rd Qu.:1.160  
##  Max.   :1.168   Max.   :1.141   Max.   :1.133   Max.   :1.214  
##       phr             mhr       
##  Min.   :0.878   Min.   :0.977  
##  1st Qu.:1.068   1st Qu.:1.013  
##  Median :1.128   Median :1.113  
##  Mean   :1.095   Mean   :1.087  
##  3rd Qu.:1.182   3rd Qu.:1.128  
##  Max.   :1.208   Max.   :1.207
##       vmrl      
##  Min.   :0.801  
##  1st Qu.:1.013  
##  Median :1.085  
##  Mean   :1.061  
##  3rd Qu.:1.128  
##  Max.   :1.193
## Highest minimum log-return: mmr
## Highest median log-return: phr
## Highest mean log-return: phr
## Highest max log-return: vhr
## cov(vmr, pmr) =  -0.001094875
## cov(vhr, phr) =  -0.0001730651

Velliv medium risk, 2011 - 2023

## 
## AIC: -27.8497 
## BIC: -25.58991 
## m: 0.0480931 
## s: 0.1198426 
## nu (df): 3.303595 
## xi: 0.03361192 
## R^2: 0.993 
## 
## An R^2 of 0.993 suggests that the fit is extremely good.
## 
## What is the risk of losing max 10 %? =< 7.4 percent
## What is the risk of losing max 25 %? =< 1.8 percent
## What is the risk of losing max 50 %? =< 0.2 percent
## What is the risk of losing max 90 %? =< 0 percent
## What is the risk of losing max 99 %? =< 0 percent
## 
## What is the chance of gaining min 10 %? >= 41 percent
## What is the chance of gaining min 25 %? >= 0 percent
## What is the chance of gaining min 50 %? >= 0 percent
## What is the chance of gaining min 90 %? >= 0 percent
## What is the chance of gaining min 99 %? >= 0 percent

QQ Plot

The qq plot looks great. Log returns for Velliv medium risk seems to be consistent with a skewed t-distribution.

Data vs fit

Let’s plot the fit and the observed returns together.

Estimated distribution

Now lets look at the CDF of the estimated distribution for each 0.1% increment between 0.5% and 99.5% for the estimated distribution:

We see that for a few observations out of a 1000, the losses are disastrous, while the upside is very dampened.

Monte Carlo

## Down-and-out simulation:
## Probability of down-and-out: 0 percent
## 
## Mean portfolio index value after 20 years: 278.454 kr.
## SD of portfolio index value after 20 years: 124.751 kr.
## Min total portfolio index value after 20 years: 0.472 kr.
## Max total portfolio index value after 20 years: 1045.731 kr.
## 
## Share of paths finishing below 100: 5.12 percent

Velliv medium risk, 2007 - 2023

Fit to skew t distribution

## 
## AIC: -34.35752 
## BIC: -31.02467 
## m: 0.05171176 
## s: 0.1149408 
## nu (df): 2.706099 
## xi: 0.5049945 
## R^2: 0.978 
## 
## An R^2 of 0.978 suggests that the fit is very good.
## 
## What is the risk of losing max 10 %? =< 5.4 percent
## What is the risk of losing max 25 %? =< 1.3 percent
## What is the risk of losing max 50 %? =< 0.2 percent
## What is the risk of losing max 90 %? =< 0 percent
## What is the risk of losing max 99 %? =< 0 percent
## 
## What is the chance of gaining min 10 %? >= 36.2 percent
## What is the chance of gaining min 25 %? >= 0.3 percent
## What is the chance of gaining min 50 %? >= 0 percent
## What is the chance of gaining min 90 %? >= 0 percent
## What is the chance of gaining min 99 %? >= 0 percent

QQ Plot

The qq plot looks good. Log returns for Velliv high risk seems to be consistent with a skewed t-distribution.

Data vs fit

Let’s plot the fit and the observed returns together.

Estimated distribution

Now lets look at the CDF of the estimated distribution for each 0.1% increment between 0.5% and 99.5% for the estimated distribution:

We see that for a few observations out of a 1000, the losses are disastrous, while the upside is very dampened. But because the disastrous loss in 2008 was followed by a large profit the following year, we see some increased upside for the top percentiles. Beware: A 1.2 return following a 0.8 return doesn’t take us back where we were before the loss. Path dependency! So if returns more or less average out, but high returns have a tendency to follow high losses, that’s bad!

Monte Carlo

## Down-and-out simulation:
## Probability of down-and-out: 0.04 percent
## 
## Mean portfolio index value after 20 years: 291.933 kr.
## SD of portfolio index value after 20 years: 120.701 kr.
## Min total portfolio index value after 20 years: 0 kr.
## Max total portfolio index value after 20 years: 2790.765 kr.
## 
## Share of paths finishing below 100: 3.38 percent

Velliv high risk, 2011 - 2023

Fit to skew t distribution

## 
## AIC: -21.42488 
## BIC: -19.16508 
## m: 0.06471454 
## s: 0.1499924 
## nu (df): 3.144355 
## xi: 0.002367034 
## R^2: 0.991 
## 
## An R^2 of 0.991 suggests that the fit is extremely good.
## 
## What is the risk of losing max 10 %? =< 8.3 percent
## What is the risk of losing max 25 %? =< 2.5 percent
## What is the risk of losing max 50 %? =< 0.4 percent
## What is the risk of losing max 90 %? =< 0 percent
## What is the risk of losing max 99 %? =< 0 percent
## 
## What is the chance of gaining min 10 %? >= 53.3 percent
## What is the chance of gaining min 25 %? >= 0 percent
## What is the chance of gaining min 50 %? >= 0 percent
## What is the chance of gaining min 90 %? >= 0 percent
## What is the chance of gaining min 99 %? >= 0 percent

QQ Plot

The qq plot looks great. Returns for Velliv medium risk seems to be consistent with a skewed t-distribution.

Data vs fit

Let’s plot the fit and the observed returns together.

Estimated distribution

Now lets look at the CDF of the estimated distribution for each 0.1% increment between 0.5% and 99.5% for the estimated distribution:

We see that for a few observations out of a 1000, the losses are disastrous, while the upside is very dampened.

Monte Carlo

## Down-and-out simulation:
## Probability of down-and-out: 0 percent
## 
## Mean portfolio index value after 20 years: 409.84 kr.
## SD of portfolio index value after 20 years: 220.813 kr.
## Min total portfolio index value after 20 years: 0.217 kr.
## Max total portfolio index value after 20 years: 2086.963 kr.
## 
## Share of paths finishing below 100: 4.16 percent

PFA medium risk, 2011 - 2023

Fit to skew t distribution

## 
## AIC: -33.22998 
## BIC: -30.97018 
## m: 0.05789224 
## s: 0.1234592 
## nu (df): 2.265273 
## xi: 0.477324 
## R^2: 0.991 
## 
## An R^2 of 0.991 suggests that the fit is extremely good.
## 
## What is the risk of losing max 10 %? =< 3.3 percent
## What is the risk of losing max 25 %? =< 0.9 percent
## What is the risk of losing max 50 %? =< 0.2 percent
## What is the risk of losing max 90 %? =< 0 percent
## What is the risk of losing max 99 %? =< 0 percent
## 
## What is the chance of gaining min 10 %? >= 32.7 percent
## What is the chance of gaining min 25 %? >= 0.1 percent
## What is the chance of gaining min 50 %? >= 0 percent
## What is the chance of gaining min 90 %? >= 0 percent
## What is the chance of gaining min 99 %? >= 0 percent

QQ Plot

The qq plot looks great. Log returns for PFA medium risk seems to be consistent with a skewed t-distribution.

##  [1] -0.091256521 -0.003731241  0.027312079  0.045808232  0.059068633
##  [6]  0.069575113  0.078454727  0.086316936  0.093536451  0.100370932
## [11]  0.107018607  0.114081432  0.127604387

Data vs fit

Let’s plot the fit and the observed returns together.

Estimated distribution

Now lets look at the CDF of the estimated distribution for each 0.1% increment between 0.5% and 99.5% for the estimated distribution:

We see that for a few observations out of a 1000, the losses are disastrous. While there is some uptick at the top percentiles, the curve basically flattens out.

Monte Carlo

## Down-and-out simulation:
## Probability of down-and-out: 0 percent
## 
## Mean portfolio index value after 20 years: 323.366 kr.
## SD of portfolio index value after 20 years: 113.866 kr.
## Min total portfolio index value after 20 years: 1.114 kr.
## Max total portfolio index value after 20 years: 3835.279 kr.
## 
## Share of paths finishing below 100: 1.76 percent

PFA high risk, 2011 - 2023

Fit to skew t distribution

## 
## AIC: -23.72565 
## BIC: -21.46585 
## m: 0.08386034 
## s: 0.1210107 
## nu (df): 3.184569 
## xi: 0.01790306 
## R^2: 0.964 
## 
## An R^2 of 0.964 suggests that the fit is very good.
## 
## What is the risk of losing max 10 %? =< 5.3 percent
## What is the risk of losing max 25 %? =< 1.4 percent
## What is the risk of losing max 50 %? =< 0.2 percent
## What is the risk of losing max 90 %? =< 0 percent
## What is the risk of losing max 99 %? =< 0 percent
## 
## What is the chance of gaining min 10 %? >= 59.6 percent
## What is the chance of gaining min 25 %? >= 0 percent
## What is the chance of gaining min 50 %? >= 0 percent
## What is the chance of gaining min 90 %? >= 0 percent
## What is the chance of gaining min 99 %? >= 0 percent

QQ Plot

The qq plot looks ok. Returns for PFA high risk seems to be consistent with a skewed t-distribution.

Data vs fit

Let’s plot the fit and the observed returns together.

Estimated distribution

Now lets look at the CDF of the estimated distribution for each 0.1% increment between 0.5% and 99.5% for the estimated distribution:

We see that for a few observations out of a 1000, the losses are disastrous, while the upside is very dampened.

Monte Carlo

## Down-and-out simulation:
## Probability of down-and-out: 0 percent
## 
## Mean portfolio index value after 20 years: 547.284 kr.
## SD of portfolio index value after 20 years: 242.954 kr.
## Min total portfolio index value after 20 years: 0.407 kr.
## Max total portfolio index value after 20 years: 1774.886 kr.
## 
## Share of paths finishing below 100: 1.26 percent

Mix medium risk, 2011 - 2023

Fit to skew t distribution

## 
## AIC: -36.9603 
## BIC: -34.7005 
## m: 0.05902873 
## s: 0.08757749 
## nu (df): 2.772621 
## xi: 0.02904471 
## R^2: 0.89 
## 
## An R^2 of 0.89 suggests that the fit is not completely random.
## 
## What is the risk of losing max 10 %? =< 3.3 percent
## What is the risk of losing max 25 %? =< 0.7 percent
## What is the risk of losing max 50 %? =< 0.1 percent
## What is the risk of losing max 90 %? =< 0 percent
## What is the risk of losing max 99 %? =< 0 percent
## 
## What is the chance of gaining min 10 %? >= 35.6 percent
## What is the chance of gaining min 25 %? >= 0 percent
## What is the chance of gaining min 50 %? >= 0 percent
## What is the chance of gaining min 90 %? >= 0 percent
## What is the chance of gaining min 99 %? >= 0 percent

QQ Plot

The fit suggests big losses for the lowest percentiles, which are not present in the data.
So the fit is actually a very cautious estimate.

Data vs fit

Let’s plot the fit and the observed returns together.

Interestingly, the fit predicts a much bigger “biggest loss” than the actual data. This is the main reason that R^2 is 0.90 and not higher.

Estimated distribution

Now lets look at the CDF of the estimated distribution for each 0.1% increment between 0.5% and 99.5% for the estimated distribution:

We see that for a few observations out of a 1000, the losses are disastrous, while the upside is very dampened.

Monte Carlo

Version a: Simulation from estimated distribution of returns of mix.

## Down-and-out simulation:
## Probability of down-and-out: 0 percent
## 
## Mean portfolio index value after 20 years: 324.124 kr.
## SD of portfolio index value after 20 years: 99.26 kr.
## Min total portfolio index value after 20 years: 3.14 kr.
## Max total portfolio index value after 20 years: 698.804 kr.
## 
## Share of paths finishing below 100: 1.2 percent

Version b: Mix of simulations from estimated distribution of returns from individual funds.

## Down-and-out simulation:
## Probability of down-and-out: 0 percent
## 
## Mean portfolio index value after 20 years: 301.082 kr.
## SD of portfolio index value after 20 years: 82.01 kr.
## Min total portfolio index value after 20 years: 64.064 kr.
## Max total portfolio index value after 20 years: 676.779 kr.
## 
## Share of paths finishing below 100: 0.32 percent

Mix high risk, 2011 - 2023

Fit to skew t distribution

## 
## AIC: -24.26084 
## BIC: -22.00104 
## m: 0.0822419 
## s: 0.07129843 
## nu (df): 89.86289 
## xi: 0.7697502 
## R^2: 0.961 
## 
## An R^2 of 0.961 suggests that the fit is very good.
## 
## What is the risk of losing max 10 %? =< 0.9 percent
## What is the risk of losing max 25 %? =< 0 percent
## What is the risk of losing max 50 %? =< 0 percent
## What is the risk of losing max 90 %? =< 0 percent
## What is the risk of losing max 99 %? =< 0 percent
## 
## What is the chance of gaining min 10 %? >= 46.1 percent
## What is the chance of gaining min 25 %? >= 1.2 percent
## What is the chance of gaining min 50 %? >= 0 percent
## What is the chance of gaining min 90 %? >= 0 percent
## What is the chance of gaining min 99 %? >= 0 percent

QQ Plot

The qq plot looks good Returns for mixed medium risk portfolios seems to be consistent with a skewed t-distribution.

Data vs fit

Let’s plot the fit and the observed returns together.

Estimated distribution

Now lets look at the CDF of the estimated distribution for each 0.1% increment between 0.5% and 99.5% for the estimated distribution:

We see that the high risk mix provides a much better upside and smaller downside.

Monte Carlo

Version a: Simulation from estimated distribution of returns of mix.

## Down-and-out simulation:
## Probability of down-and-out: 0 percent
## 
## Mean portfolio index value after 20 years: 498.255 kr.
## SD of portfolio index value after 20 years: 157.588 kr.
## Min total portfolio index value after 20 years: 133.651 kr.
## Max total portfolio index value after 20 years: 1462.848 kr.
## 
## Share of paths finishing below 100: 0 percent

Version b: Mix of simulations from estimated distribution of returns from individual funds.

## Down-and-out simulation:
## Probability of down-and-out: 0 percent
## 
## Mean portfolio index value after 20 years: 479.55 kr.
## SD of portfolio index value after 20 years: 163.678 kr.
## Min total portfolio index value after 20 years: 62.493 kr.
## Max total portfolio index value after 20 years: 1273.474 kr.
## 
## Share of paths finishing below 100: 0.1 percent

Compare pension plans

Risk of max loss of x percent for a single period (year).
x values are row names.

Velliv_medium Velliv_medium_long Velliv_high PFA_medium PFA_high mix_medium mix_high
0 21.3 18.2 19.9 12.2 14.3 12.7 13.0
5 12.5 9.6 12.8 6.0 8.6 6.2 4.2
10 7.4 5.4 8.3 3.3 5.3 3.3 0.9
25 1.8 1.3 2.5 0.9 1.4 0.7 0.0
50 0.2 0.2 0.4 0.2 0.2 0.1 0.0
90 0.0 0.0 0.0 0.0 0.0 0.0 0.0
99 0.0 0.0 0.0 0.0 0.0 0.0 0.0

Chance of min gains of x percent for a single period (year).
x values are row names.

Velliv_medium Velliv_medium_long Velliv_high PFA_medium PFA_high mix_medium mix_high
0 78.7 81.8 80.1 87.8 85.7 87.3 87.0
5 63.8 64.9 69.2 71.5 75.8 71.4 69.9
10 41.0 36.2 53.3 32.7 59.6 35.6 46.1
25 0.0 0.3 0.0 0.1 0.0 0.0 1.2
50 0.0 0.0 0.0 0.0 0.0 0.0 0.0
100 0.0 0.0 0.0 0.0 0.0 0.0 0.0

MC risk percentiles: Risk of loss from first to last period.
_a is simulation from estimated distribution of returns of mix.
_b is mix of simulations from estimated distribution of returns from individual funds.

_m is medium.
_h is high.

Velliv_m Velliv_m_long Velliv_h PFA_m PFA_h mix_m_a mix_h_a mix_m_b mix_h_b
0 5.12 3.38 4.16 1.76 1.26 1.20 0 0.32 0.10
5 4.52 2.98 3.84 1.54 1.10 1.04 0 0.24 0.06
10 3.76 2.42 3.48 1.46 1.06 0.96 0 0.18 0.06
25 2.48 1.54 2.52 0.98 0.72 0.70 0 0.08 0.04
50 0.88 0.72 1.16 0.52 0.38 0.40 0 0.00 0.00
90 0.06 0.12 0.20 0.10 0.08 0.02 0 0.00 0.00
99 0.04 0.04 0.04 0.00 0.04 0.00 0 0.00 0.00

MC gains percentiles: Chance of gains from first to last period.
_a is simulation from estimated distribution of returns of mix.
_b is mix of simulations from estimated distribution of returns from individual funds.

Velliv_m Velliv_m_long Velliv_h PFA_m PFA_h mix_m_a mix_h_a mix_m_b mix_h_b
0 94.88 96.62 95.84 98.24 98.74 98.80 100.00 99.68 99.90
5 94.20 96.06 95.50 98.06 98.62 98.66 100.00 99.60 99.88
10 93.18 95.84 95.04 97.88 98.46 98.46 100.00 99.48 99.82
25 90.92 94.04 93.52 96.98 98.24 98.04 100.00 99.04 99.64
50 85.08 90.48 90.80 94.88 97.24 96.18 99.88 97.62 99.28
100 71.16 78.74 83.04 88.50 94.64 89.84 99.34 89.70 97.40
200 39.22 43.14 65.08 59.88 84.98 58.64 92.50 48.42 87.12
300 15.58 16.88 45.62 22.60 70.72 21.98 71.48 11.38 66.20
400 5.26 4.30 30.20 4.04 53.34 3.88 44.12 1.24 41.32
500 1.50 1.14 18.38 0.44 37.90 0.30 23.12 0.08 21.70
1000 0.00 0.02 0.58 0.02 2.24 0.00 0.38 0.00 0.14

Summary statistics

Fit summary

Summary for fit of log returns to an F-S skew standardized Student-t distribution.
m is the location parameter.
s is the scale parameter.
nu is the estimated degrees of freedom, or shape parameter.
xi is the estimated skewness parameter.

Velliv_medium Velliv_medium_long Velliv_high PFA_medium PFA_high mix_medium mix_high
m 0.048 0.052 0.065 0.058 0.084 0.059 0.082
s 0.120 0.115 0.150 0.123 0.121 0.088 0.071
nu 3.304 2.706 3.144 2.265 3.185 2.773 89.863
xi 0.034 0.505 0.002 0.477 0.018 0.029 0.770
R-squared 0.993 0.978 0.991 0.991 0.964 0.890 0.961

Monte Carlo simulations summary

Monte Carlo simulations of portfolio index values (currency values).
Statistics are given for the final state of all paths.
Probability of down-and_out is calculated as the share of paths that reach 0 at some point. All subsequent values for a path are set to 0, if the path reaches at any point.
0 is defined as any value below a threshold.
losing_prob_pct is the probability of losing money. This is calculated as the share of paths finishing below index 100.

## Number of paths: 5000
Velliv_m Velliv_m_long Velliv_h PFA_m PFA_h mix_m_a mix_m_b mix_h_a mix_h_b
mc_m 278.454 291.933 409.840 323.366 547.284 324.124 301.082 498.255 479.550
mc_s 124.751 120.701 220.813 113.866 242.954 99.260 82.010 157.588 163.678
mc_min 0.472 0.000 0.217 1.114 0.407 3.140 64.064 133.651 62.493
mc_max 1045.731 2790.765 2086.963 3835.279 1774.886 698.804 676.779 1462.848 1273.474
dao_prob_pct 0.000 0.040 0.000 0.000 0.000 0.000 0.000 0.000 0.000
losing_prob_pct 5.120 3.380 4.160 1.760 1.260 1.200 0.320 0.000 0.100
## Highest mean    : PFA_h
## Lowest sd       : mix_m_b
## Highest min     : mix_h_a
## Highest max     : PFA_m
## Lowest dao prob : Velliv_m Velliv_h PFA_m PFA_h mix_m_a mix_m_b mix_h_a mix_h_b
## Lowest loss prob: mix_h_a

Appendix

Average of returns vs returns of average

Definition: R = 1+r

## Let x_0 be 100.
## Let y_0 be 200.
## So the initial value of the pf is 300 .
## Let R_x be 0.5.
## Let R_y be 1.5.

Then,

## x_1 is R_x * x_0 = 50.
## y_1 is R_y * y_0 = 300.

Average of returns:

## 0.5 * (R_x + R_y) = 1

So here the value of the pf at t=1 should be unchanged from t=0:

## (x_0 + y_0) * 0.5 * (R_x + R_y) = 300

But this is clearly not the case:

## 0.5 * (x_1 + y_1) = 0.5 * (R_x * x_0 + R_y * y_0) = 175

Therefore we should take returns of average, not average of returns!

Let’s take the average of log returns instead:

## 0.5 * (log(R_x) + log(R_y)) = -0.143841

We now get:

## (x_0 + y_0) * exp(0.5 * (log(Rx) + log(Ry))) = 259.8076

So taking the average of log returns doesn’t work either.

Simulation of mix vs mix of simulations

Test if a simulation of a mix (average) of two returns series has the same distribution as a mix of two simulated returns series.

We are adding annual returns rather than multiplying, so imagine that we are simulating log returns.

## m(data_x): -0.08587544 
## s(data_x): 0.3849157 
## m(data_y): 10.16474 
## s(data_y): 2.865587 
## 
## m(data_x + data_y): 5.03943 
## s(data_x + data_y): 1.397076

m and s of final state of all paths.
_a is mix of simulated returns.
_b is simulated mixed returns.

m_a m_b s_a s_b
100.791 100.755 6.728 6.210
100.894 100.715 6.500 6.291
100.876 100.517 6.441 6.417
101.060 100.761 6.404 6.176
100.763 100.925 6.267 5.989
100.675 101.274 6.546 5.998
100.915 100.640 6.466 6.255
100.664 100.595 6.554 6.216
101.150 100.909 6.586 6.156
100.734 100.877 6.361 6.010
##       m_a             m_b             s_a             s_b       
##  Min.   :100.7   Min.   :100.5   Min.   :6.267   Min.   :5.989  
##  1st Qu.:100.7   1st Qu.:100.7   1st Qu.:6.413   1st Qu.:6.046  
##  Median :100.8   Median :100.8   Median :6.483   Median :6.193  
##  Mean   :100.9   Mean   :100.8   Mean   :6.485   Mean   :6.172  
##  3rd Qu.:100.9   3rd Qu.:100.9   3rd Qu.:6.552   3rd Qu.:6.246  
##  Max.   :101.2   Max.   :101.3   Max.   :6.728   Max.   :6.417

_a and _b are very close to equal.
We attribute the differences to differences in estimating the distributions in version a and b.

The final state is independent of the order of the preceding steps:

So does the order of the steps in the two processes matter, when mixing simulated returns?

The order of steps in the individual paths do not matter, because the mix of simulated paths is a sum of a sum, so the order of terms doesn’t affect the sum. If there is variation it is because the sets preceding steps are not the same. For instance, the steps between step 1 and 60 in the plot above are not the same for the two lines.

Recall, \[\text{Var}(aX+bY) = a^2 \text{Var}(X) + b^2 \text{Var}(Y) + 2ab \text{Cov}(a, b)\]

## [1] 0.005355618
## [1] 0.005355618

Our distribution estimate is based on 13 observations. Is that enough for a robust estimate? What if we suddenly hit a year like 2008? How would that affect our estimate?
Let’s try to include the Velliv data from 2007-2010.
We do this by sampling 13 observations from vmrl.

##        m                 s          
##  Min.   :0.05922   Min.   :0.04444  
##  1st Qu.:0.06576   1st Qu.:0.06218  
##  Median :0.06906   Median :0.06761  
##  Mean   :0.07029   Mean   :0.06849  
##  3rd Qu.:0.07396   3rd Qu.:0.07422  
##  Max.   :0.08599   Max.   :0.09295